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ABSTRACT 
The United States Navy's (USN) numerical front~location parameter, 
related to relative maxima in thermal-wind shear, as developed by Renard 
and Clarke, is investigated in its relation to the data-processing system 
used at the USN Fleet Numerical Weather Facility (FNWF), Monterey, 
California and to three front-associated wind parameters: a) geostrophic 


wind component normal to isentropes, b) shear of geostrophic wind 
& 


Fy 


component tangent to isentropes, and c) m frontogenetical parameters, 


‘ 


eg. 
e 


The study, which generalizes results to hyperbaroclinic zones and their 
boundaries, is carried out for the Northern Hemisphere at 1000, 850, 700, 
500, and 300 mb for 00Z and 12Z, 20 August 1965, in the framework of 
FNWF's numerical analysis system. 

Results indicate some major differences between the sytioptic-scale 
baroclinic zones analyzed from radiosonde~observed temperatures, and 
FNWF's numerically processed mandatory~level virtual temperature data, 
at 1000 mb, over ocean areas, and at latitudes north of 55 degrees 
latitude. The normal geostrophic wind component was found to be useful 
in designating the type of front (cold, warm, stationary). Frontal shear 
proved to be mostly cyclonic at low levels (1000, 850 mb) becoming 
increasingly anticyclonic at higher levels (700, 500, 300 mb). The 
advective frontogenesis parameter appeared to be of limited value only. 
Where appropriate the numerical products are compared to the National 


Meteorological Center's frontal analyses. 
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1. Introduction 

The need for an objectively specified front locator has 
prompted interested individuals at the United States Navy Fleet 
Numerical Weather Facility (FNWF) and United States Naval 
Postgraduate School (NPGS) to attempt a solution by numerical 
methods. (3] In their search for such a parameter Renard and 
Clarke tried various schemes. The one giving best results involves 
the second derivative of potential temperature. GG@, the name 
given to the frontal parameter, is related to the shear vorticity of 


the thermal wind and may be symbolized as_ [I] 
~Vivel*VO _ _ oO ~<a 
= Vive tle = 9%, (1) 


Vel 
GG@, as computed at FNWF, is based upon potential 





temperature (@) determined from a unique radiosonde-data processing 
system which employs heights of seven mandatory pressure levels 
(Zp), and sea-level pressure. [s] A static stability parameter 
(SS), expressed as a function of Zp, is then derived for five 

layers (1000-775 mb; 775-600 mb; 600-450 mb; 450-300 mb; 

300-200 mb). The SS is assumed constant in each layer, but is 
allowed to vary along a constant-pressure surface. Finally, 

virtual temperatures are computed at 3969 Northern Hemisphere 

grid points using the 1000-mb height, 1000/500 mb thickness 

and SS. Thus, the derived 6 is a virtual potential temperature 


Sona 


with the resulting fields hydrostatically consistent, statically 


Jk 


stable, and horizontally compatible. This particular @ (or T) is 
hereafter referred to as a processed @ (or T), 

It is to be noted that although GG@ is used operationally to 
locate fronts numerically, the parameter is, in reality, a 
generalized baroclinic-zone locator. Figure 4 shows a typical 
FNWF analysis of GG@, in this case the 850-mb level. Ridges 
and troughs in the GG@ field represent warm and cold boundaries 
of hyperbaroclinic zones, respectively, and as such are axes of 
maximum shear in the thermal wind. Since GG@ has been 
observed to have a direct relation to VV 6, the intensity of GG@ 
at the baroclinic-zone boundaries gives an indication of the 
maximum strength of the baeanitic zone. 

The warm-air boundaries of the hyperbaroclinic zones 
represent the numerical counterparts of conventional fronts. 
Figure 14 shows a typical relation between FNWF's numerically 
determined fronts (xxxx) and the USWB's hand-analyzed fronts 
(- - -). As may be noted, there is often a close comparison, 
espegially in dense data areas. In other areas the two fronts 
are quite incompatible in orientation and/or location. Recent 
studies by Clarke and Renard’! ‘and Hamrick have shown a need 
for improvement of the numerical front so as to be more in tune 


with hand-analyzed fronts. [1, 4] 
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2. Objectives and Outline of Study 

This thesis study concentrated on two as yet unexplored 
areas of numerical-front analysis; 1) the effect of numerical 
processing on location and intensity of numerical fronts; and 
2) the relation of pertinent wind parameters to the frontal analysis. 

A comparison of the processed virtual temperature with 
reported radiosonde temperature was conducted by the author for 
the Northern Hemisphere at five levels: 1000 mb, 850 mb, 700 mb, 
500 mb, and 300 mb. A similar comparison was also made for 
the GG@ fields derived from the two temperature fields. 

Fronts as used operationally by synoptic meteorologists are 
associated with unique wind phenomena. Presently, FNWF's 
numerical~front model does not incorporate any wind parameters. 
In this paper three such wind parameters, as determined 
numerically, are investigated in association with both the 
computerized and manually~analyzed fronts and baroclinic zones. 

The first of these wind parameters considered was the 
component of the geostrophic wind normal to the isentropes (VNF) 
on an isobaric surface, where the isentropes, of necessity, are 
regarded as a best approximation to the orientation of the fronts. 


This is expressed as: 
~\7-¥9 ~i7.y 
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Since the motion of a front is generally taken as the component 
of wind normal to it, the movement of numerical fronts should be 
statistically related to VNF. 

second, the frontogenical properties of a numerically 


determined front may be indicated by: [ 6] 
FG =-V(%" hal) VO =~ 0(G-2)-Ve 
=-V(VAF)*V@ 


Since GG@ is well correlated with baroclinicity, time changes 
in the magnitude of GG®@ may be used as a measure of frontogenesis 
in order to verify the applicability of equation (3). 

The third wind parameter betterueuse was shear of the 
geostrophic wind component parallel to the isentropes on an 


isobaric surface: 
_ 7 . + | 2) eee 49) 
SHVF vv K x ve | ve 
— A N iD 
sy \G° K x A, ) “(e 


Horizontal wind shear is a parameter often used to locate 


(4) 


conventional fronts, and so the relation of shear to both the 
baroclinic zone and its warm-air boundary is considered here. 
The study of these various wind parameters was carried out 
for all five of the previously mentioned levels, for both the FNWF 
processed temperatures and the observed temperatures at two 


radiosonde reporting times. 
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3. Data and Charts 

Data and chart material for the whole Northem Hemisphere, 
collected at 0000Z and 1200Z 20 August 1965, were utilized for 
this study. The selection of these times was purely arbitrary in 
that no unseasonable weather conditions influenced their choice. 
For certain portions of the study, only the dense~data areas of 
the North American Continent and North Atlantic Ocean were 
utilized. The basic chart used is the 1:30,000,000 polar stereo- 
graphic Northern Hemisphere projection with data at 63 x 63 grid 
points, each 38]lkm apart at 60 N. 

The chart material consisted of 6, V6, GG@, VNF, FG, and 
SHVF obtained from both the FNWF processed virtual temperature 
and observed temperature data. The latter were obtained from 
FNWF listings of mandatory-level temperature data. These 
observed temperatures were extracted from the data printouts, 
put on punch cards and numerically analyzed at FNWF, using the 
processed temperature fields as first guess. 

Rather extensive procedures were adopted to uncover errors 
in transcribing data from listings and punch cards. The fields 
of analyzed temperatures and differences in processed and 
observed temperatures were carefully screened for indication of 
errors in the observed temperature fields. 

Numerical differentiation of wind and temperature parameters 
was accomplished by a finite differencing technique employing the 


quartic interpolation polynomial in the form 
lies, 


A\x,y( ) At 2e=0 
-izal hy r§l Me “0):.)) -( Jee } (5) 


After the 12Z 20 August 1965 observed temperature charts had 
been analyzed for frontal phenomena, it was discovered that an 
improper guess field had been used. A reanalysis of these 
temperatures using a corrected guess field, showed insignificant 
changes from the initial temperatures for grid points north of 15 N, 
so the 12Z chart material as originally run, was: accepted. 

4. Temperature: FNWF Processed vs. Observed 

To determine the difference between the FNWF processed 
(virtual) temperature and the observed temperature, the latter 
was subtracted from the former at each grid point. Since the 
FNWF temperatures are virtual, while the observed are not, 
differences of 0.0 C to +3.0 C are to be expected. 

The resulting temperature differences were categorized by 
level, locale, and magnitude. The stratification by location 
involved a latitudinal grouping (South, 20 N r 35 N; Central, 

35 N - 55 N; and North, 55 N - 90 N), and a longitudinal grouping 
(North America, 70 W - 130 W; Pacific, 130 W - 140 E; Asia, 

140 E - 60 E; Europe, 60 E -— 10 W; Atlantic, 10 W - 70 W). 

Figure 1 shows the limits of these various areas. The number 

of grid points within each area is given alongside the title of 


each area. 
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Temperature differences for 00Z 20 August 1965 are 
tabularized by per cent of grid points (see Table 1). It can be 
noted from the table that the South area has the greatest per cent 
of grid points in the 0.0 to +2.9 range and the least in the -0.9 
to 0.0 column at all levels. Except for 500 mb, the reverse is 
true of the North area. These figures merely reflect the normal 
meridional gradient of moisture and give reason for the over-all 
larger V0 and GG@ using the FNWF vice the observed temperatures. 
Differences =+3.0C and £ -0.9 C are not explainable in terms 
of moisture nor may they be regarded as negligible “noise” 
differences. Such cases amount to an average of 11, 12, and 20% 
in the South, Central, and North zones,respectively ,for all levels 
except 1000 mb, while 1000 mb shows 37, 55, and 44% of cases 
in the three zones, respectively. 1000 mb, of course, iS more 
often than not a fictitious atmospheric level and this may account 
forthe discrepancies, while such things as mislocated data, 
transmission errors, processing errors, and so forth, account for 
the remainder, 

When the longitudinal division of area is employed, the high 
moisture content of the air above water surfaces is evidenced by 
the higher percentages over the Pacific and Atlantic in the 0.0 
to +2.9 range, especially at 850 and 1000 mb (Table 2). As 
already indicated above, 1000 mb shows the maximum differences, 


while the Asiatic area is noted to have the greatest percentage of 


large temperature differences at all levels. 


vy 


Thus, collectively, the two tables suggest that GG@p and 
GG6o fields are most likely to differ at the 1000*mb level, at 
latitudes poleward of 55 N, and over the Asiatic area 60 E- 140 E. 

Ferrentino, in a climatological study, made a similar 
comparison of temperature differences for the same observation 
time at the 200, 500, and 850-mb levels. (3] He initially 
converted the observed temperatures to virtual temperatures and 
calculated differences from the FNWF-processed temperatures 
at reporting stations only. His findings indicated that the 500é¢mb 
level was the best of the three levels tested, with 80 per cent of 
the stations having a difference within plus of minus 1.2C. This 
result its in good agreement with Tables 1 and 2, which show the 
least percentage of unexplainable large temperature differences 
at 500 and 700 mb. 

9. GG6 Analyses 


Criteria Used in the Analysis of 
GG6 Isolines, Ridges, and Troughs 


Values of GG@ were computed and printed out on the FNWF 
63 x 63 square grid (as were all other parameters being considered) 
in units of 1072C/(100km)¢. Both positive and negative GG@ 
areas are indicated (Fig. 2-11) by isolines drawn on Northern 
Hemisphere 1 : 30,000,000 polar stereographic projections at 
intervals of 10 units starting with the 5 line. GG@® values are 
generally less than 35 units. Closed contours were not drawn 


around single grid point values =|5|units of GG@. 


18 


The main requirement for the placement of fronts (i.e., GG@ 
ridges or warm-air boundaries of baroclinic zones) on the 
contoured GG@ fields was that the front be associated with a 
baroclinic zone that had a well defined rear (cold~air) boundary. 
Thus, in order for a ridge in the GG6@ field to be indicated, it 
must have an associated trough, with both ridge and trough 
meeting certain necessary requirements, as follows. Ridges and 
troughs were only drawn within an area of GG@ where magnitudes 
exceeded five units. An exception to this rule occurred when two 
distinct ridge lines, aligned end to end, were separated by a 
small distances of léss than two imeshilengthsihavirlg valués between 
zero and +5 units. This exception was not allowed if the trough 
line was also discontinuous in the same fashion. This rule was 
also made applicable to trough lines if they were Separated by a 
similarly small distance of negative values. Areas of positive 
and negative GG@ that qualified, but were circular in shape, such 
that distinct ridges or troughs could not be discerned, were also 
not analyzed for these entities. No distance limitations between 
ridge and trough were used, but ridges and troughs were required 
to be approximately parallel. 

The reason for the preceding limiting conditions was to 
specify an objective means of locating only the synoptic~-scale 
fronts. It can also be seen that if all GG6@ ridges, without 


further refinement, were classified as fronts, the number and 
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mileage of numerically-determined fronts would far exceed that 
for operationally produced hand=-analyzed fronts, resulting in a 
situation not entirely desirable for this study. 

Comparison of GG@p, GGeo, and USWB Fronts at 00Z 20 Aug 1965 

A descriptive comparison of the two GGé@ fields, derived 
from FNWF processed (GG@p) and observed (GG@o) temperatures, 
for 00Z at five mandatory levels, is given below. 

(1) 1000 mb (Eligisg2).3) 

The fronts indicated by the two GG@ fields are basically in 
very poor agreement, a situation somewhat expected from the 
results shown in Tables 1 and 2. The GG@fields show a very 
limited number of closed centers over ocean areas, while the GG@p 
field has centers over ocean areas, but decreased intensity 
compared to land areas. The fronts along the western part of 
North America, for which there is some correlation between the 
two charts, &arei probably a result of diurnal heating and the trade 
wind inversion, a phenomena peculiar to the West Coast, 
especially during summer. 

(2) 850 mb (Figs. 4, 5) 

The locations and intensities are in good agreement, especially 
over eastern United States. Relation of the two fields is 
particularly poor poleward of 55 N in the North American Arctic. 
Intensities of GG®remains less than GG@p over oceans. Inversion 
and diurnal heating effects appear to continue domination of the 


North American west coast. 
20 


(3) 700 mb (Figs. 6!).7) !) 

This level shows more frontal mileage than lower levels, 
mostly due to increases over ocean areas. Agreement of fronts 
has improved over 850 mb. The United States west coast 
baroclinic zone has changed orientation. 700 mb is normally 
above the trade wind inversion. 

(4) 500 mb (Figs. 8, 9) 

Overall intensity on both charts is large with GG@p values 
greater than GG@o. Agreement of position and number of fronts 
is better than at lower levels. Whereas lower levels featured 
short~segmented fronts, there are now extensive frontal bands. 
Low latitude frontal bands are now quite evident. 

(S) 300 mb (Pigs .20, 11) 

A marked increase of intensity compared to lower levels is 
noted for each field; also fronts are more continuous. The greatest 
variation in front location from Figures 10 and 11 is shown over 
China. The low-latitude baroclinic zone system is well developed. 

The position of the fronts at the five levels,for processed and 
observed temperatures,are given in Figures 12 and 13, respectively. 
For processed temperatures, the frontal contours that lie over the 
eastern United States, Canada, northern Atlantic, northwestern 
Europe, and Be nizal Asia are internally consistent. Slopes range 
from 50 : 1 in the cold front sector to 150 : 1 in the warm front 
sector. Such correspondence is not as evident in the observed 
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frontal contours. In areas other than those mentioned above, a 
display of vertically consistent slopes on both charts is not so 
evident. Numerous occasions of intersecting frontal contours 

and fronts sloping toward warm air with increasing height are noted, 
The two figures make obvious the poor fit of 1000;mb fronts with 
those at other levels. 

When GG®@oiand GG@p fronts are compared with USWB fronts 
for 850 mb (Fig. 14), a closer comparison between GG@o and GG@p 
fronts exists than between USWB fronts and either of thérGGé fronts 
Incompatabilities occur for USWB occlusions and east ocean~area 
fronts, and GG@ fronts over wester North America and Africa. 

Comparison of Frontal Positions at 00Z and 12Z 20 Aug 1965 

Figure 15 shows the position of GG@p fronts and USWB fronts 
for 00Z and 12Z. In areas where the two may be associated, such 
as easter North America, western Europe, andthe west Pacific 
similar 12-hour frontal movements are noted. 

Summary 

In summarizing the differences that occur between the two GG@ 
fields, the weak GG6o baroclinic zones over water for the lower 
levels were most obvious. The GG@p zones are relatively 
stronger due to processed temperatures being virtual temperatures. 
In effect, when crossing from the rear to the front of a transition 
zone, that is from cold dry air to warm moist air, the temperature 


range for observed temperatures is less than that for processed 


if 


temperatures. This reduction of temperature range directly 
decreases the intensity of GG@6o. At upper levels, the air contains 
less moisture, and consequently this variation between GG8o 

and GG@p does not occur. 

Fronts found at 700 mb and at higher levels were indicated by 
both GG@ fields. The existence of upper-level fronts has been 
explained as intrusions of stratospheric air into the troposphere 
through breaks in the tropopause. [7] These intrusions were 
observed in polar air. The upper-level GG@ fronts,which ogcur 
in tropical air, are possibly associated with the tropopause break 
surrounding the sub-tropical jet stream. 

6. VNF Analyses 
Ee tON 

In applying equation (2) it was assumed that isentropes are 
aligned with GG@ ridges. A check at all levels, for both processed 
and observed temperatures, showed this to be a very good 
approximation. With reference to Figure 16, showing a cold=front 
situation at 700 mb, the following explains schematically the 
positive and negative zones in the VNF field. Here, since the 
geostrophic component, VNF, is directed along Ne , VNF is 
positive. For warm front case, VNF and Ne are oppositely 
directed and hence VNF is negative. 

The relation of actual frontal movement to the component of 
geostrophic wind normal to isentropes, formed from both processed 
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temperatures (VNFp) and observed temperatures (VNFo), for the 
period 0000Z to 1200Z, was studied. Figure 17 shows a typical 
VNFp pattern over North America and the North Atlantic at 850 mb 
for 0000Z, with contours at 0, +2.5, +5.0, and +10.0 meters/sec. 
Warm and cold air boundaries of baroclinic zones are also 
indicated, so that their relation with the VNF field may be noted. 
Fronts of any significant length are seen to be crossing repeatedly 
through positive and negative areas of VNFp, which, in general, 
are cellular in shape. 


Comparison of VNFp to GG@p Baroclinic Zones 
at 00Z and 12Z 20 August 1965 


The relation of VNFp to the processed=-temperature GGé@ fronts, 
at the five mandatory levels for the 12-hour period, will be discussed 
first. Since 850 mb offers a good description of the existing large- 
scale conditions near the surface, without undue modification 
resulting from surface influences, it will be discussed in greater 
detail. Variations from 850 mb phenomena will be emphasized in 
the discussion of the other levels. 

(1) 850 mb (rigem 7, Foe 19) 

Figures 17 and 18 show the VNFp pattern and the fronts for 
times 00Z and 12Z, respectively. An alternation of positive and 
negative VNFp centers along the front occur at both times. As 
expected from the VNF pattern along the baroclinic zone, the 
amplitude of the frontal wave increases during the period (see 
Fig. 19). 
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A traveling frontal wave is indicated by the movement of 
negative and positive VNFp centers during the period. To wit, 
the positive area of VNFp, having its greatest intensity over 
James Bay in Canada at 00Z, moves south southeastward to within 
the baroclinic zone over northern New York state by 12Z; the 
negative center,just off Nova Scotia and Newfoundland at 00Z, 
moves northeastward to just east of Newfoundland by 12Z. This 
movement, which averaged close to 40 knots, and movement of 
other pertinent VNFp centers, is shown by history positions on 
the 12Z chart. 

All VNFp positive and negative areas having at least one 
closed contour at 00Z were found to have counterparts at 12Z, 
except for two centers located in the extreme northeastern comer 
of the area, which most likely moved out of our area of study. 
New centers can sometimes be traced to previous centers that have 
split. 

The relation of VNFp to frontal movement was checked by 
selecting points within the hyperbaroclinic zone and then 
calculating the average VNFp and actual movement experienced 
by these identifiable points during the 12-hour period. Of nine 
uniformly distributed points, eight moved in the direction indicated 
by the 12-hour average of VNFp. The magnitude of VNFp was also 


a good indicator of the amount of movement. 


Zo 


(2) 1000 mb (Rigisa 207221) 

The VNFp field for this level has the characteristic alternating 
VNFp centers along the front. The vertical consistency between 
1000 mb and 850 mb is quite good at both 00Z and 12Z, except 
for a small area just off Nova Scotia. The intensity of the closed 
centers at 1000 mb is slightly below the intensity of the 850 mb 
centers. Time consistency at the former is less favorable than 
that of 850 mb. The relation of VNFp to movement of the baroclinic 
zone showed that of six baroclinic-zone points checked, only 
three showed direction of movement indicated by VNFp. 

(3) 700 mb (tio Figs. shown) 

Counterparts of all major centers at 850 mb are observed 
at 700 mb, with little change in intensity or shape for both 
observation times. Time consistency for this level is good in 
spite of an increased amount of splitting and combining of centers 
during the 12-hour period. Negative centers tended to be more 
intense than positive centers. Seven out of ten points in the 
transition zone showed agreement of VNFp to baroclinic zone 
movement. 

(4) 900 mb (Figs. 22,923) 

The time consistency for this level has deteriorated from that 
of lower levels. In conjunction, the degree of space continuity 
with lower levels has fallen off. The intensity values of centers 
has gone unchanged from 700 mb. VNFp was an indicator of 


direction of movement at only five of eight baroclinic~zone:points 


investigated. 
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(5) 300 mb (no Figs. shown) 

A marked increase of intensity of VNFp centers relative to 
lower levels is the result of a sharp increase in the Vg, and not 
a more normal alignment of Vg to to the isentropes. An increase 
of centers with decreased sizes was also noted. The tracing of 
VNFp centers from one observation time to the next is more difficult, 
and only partial space continuity exists. VNFp is not a good 
indicator of front movement due to probable non-linear fluctuations 
of VNFp within the observation period. 


Comparison of VNFo to GG6o Baroclinic Zones 
at 00Z and 12 Z 20 August 1966 


The comparison of VNFo with the observed temperature GG@ 
front movement is presented for two levels only, 850 mb and 500 mb. 

(1) 850 mb (Figs. 24, 25, 26, 27) 

Figures 24 and 25 show the VNFo field in relation to the warm- 
air baroclinic-zone boundaries for the two observation times. The 
alternation of VNFo centers along the front is similar to that 
detected for the processed temperature charts (Figs. 17, 18, 19). 
Since the same Vg field is used for both VNFs (processed and 
observed) for a given time, differences between the two must result 
from differences in the two VY @ fields. Figure.26 shows a 
comparison of the two VNF fields for 00Z. The only major differences 
are those near Greenland and 6ver:the extreme northea'st Atlantic. 

As shown by Figure 27, there is an outstanding time consistency 


of VNFo centers for the time period. Positive centers generally 


ad 


show a south and east movement, while negative centers show a 
north and east movement, similar to that observed for processed 
temperature VNF centers. However, the expected increase in 
amplitude of the 122 fronts does not materialize, although most 
of the North American-North Atlantic frontal systems do move in 
accordance with the VNFo fields. 

(2) 500 mb (no Figs. shown) 

The time continuity of the VNFo centers has fallen off 
considerably when compared to the ddwerlevels. Space continuity 
from the 850-mb level is quite limited. Only along small sections 
of the front can VNFo be considered an accurate indicator of front 
movement. 

Comparison of VNF Fields with USWB Fronts 

When VNF fields derived from processed and observed 
temperatures were compared with USWB fronts, only the 850-mb 
level was considered. This limitation was due to the fact that 
the USWB operationally produces only 850-mb and surface analyses 
of fronts. 

The VNFp centers were excellent indicators of the instantaneous 
movement of USWB fronts for all portions except for the extension 
into Texas (Figs. 28, 29). The USWB front at 00Z has a wave 
structure which increased in amplitude; eastward progression of 
this wave was well indicated by the location and movement of the 
VNFp centers. This is shown in fine detail along the front in the 


North Atlantic. 
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When the VNFo field was compared to the USWB tronts for 
the 850-mb level, a high degree of relationship prevailed (Figs. 30, 
31); however, it was not as profound as that obtained using the 
processed~-temperature VNFs. 

summery 

In summary, it was noted that when a front came under the 
influence of an alternating VNF field, the front tended to develop 
a progressive wave-like structure. This phenomena, which is 
known to exist in nature, was most evident at 850 mb, using 
processed temperature. The VNF fields for both processed and 
observed temperatures had time continuity at each level and space 
continuity between levels, except at 300 mb. This upper level 
also had comparatively smaller centers along with an increase in 
intensity. When the 850-mb VNF was applied to the USWB fronts 
for this level, a high degree of correlation with the front movement 
existed, especially when using processed temperatures. 
7. FGAnalyses 

The expression used for frontogenesis (FG) here (equation 3) 
is an advective one, which in effect, assumes conservation of 
temperature by individual air particles. Values of FG were 
computed numerically at grid points and expressed in units of 0.781 
x 107°C/(100km)/sec. Figure 32 shows the results obtained at 
850 mb using processed temperatures, The larger values of FG are 


situated within the hyperbaroclinic zone, and rightly so, as FGis 
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proportional to V/@. Intense centers of FG, at all levels, showed 
good time continuity. This was true of space continuity also, 
which was best between 850 and 700 mb. 

It is to be noted that the FG values are instantaneous 
quantities, however time-integrated values are of interest for 
comparison to changes in the GG@ field. In order to assess the 
relation of FG to changes in GG@ along the numerical front, it 
was assumed that if a grid point was experiencing positive FG 
at both the beginning and the end of a 12-hour period, it 
experienced positive FG throughout. The same is assumed of grid 
points having negative values (i.e., frontolysis). Ifa grid point 
experienced a negative (positive) value initially and a positive 
(negative) value at the end of the period, it was not considered 
in the evaluation. Figure 33 is a 1000-mb depiction of North 
Americanand North Atlantic areas having oonstarit-algebraiaosign of FG 
during the 12-hour period, 00Z to 12Z 20 August 1965. It can be 
seen that about two thirds of the area showed no change in 
algebraic sign of FG. This percentage, to within 10 per cent, 
applies to all other levels, except 300 mb, which had a markedly 
decreased percentage. 

The evaluation is summarized in Table 3, which shows the 
per cent of FNWF grid points at which time variations in Y @ agree 
and disagree with the FG calculations. Frontogenesis (frontolysis) 
Occurs when Y@ increases (decreases) during the period. 300 mb 
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is not included due to the numerous points showing changes in 
sign of FG during the period, thereby giving few grid points for 
verification. 

The results given in Table 3 imply that the computed positive 
FG is not a good indication of frontogenesis, although negative 
FG (i.e., frontolysis) is more promising, especially at the lower 
free-atmosphere levels of 850 and 700 mb. An explanation of this 
may be that a certain amount of frontogenesis is needed to maintain 
a mature baroclinic zone against such non-~advective frontolytical 
processes as radiation, conduction, convection, and mixing, 
whereas frontolysis is associated with a definite decrease in the 
intensity of YW@as the non~advective processes are frontolytical 
in Character. With regard to FG indicating formation as well as 
intensification, the FG pattern at 500 mb for 00Z was indeed a 
good indicator of baroclinic-zone formation by 12Z in the North 
Atlantic. But in other areas the relation was not so good. 

The better results at 850 and 700 mb give evidence that 
advective frontogenical processes are the prime factors at lower 
levels as noted by Sanders. [9J But then why doesn't 1000 mb 
have the best verification? This may result from the geostrophic 
wind being a poor representation of the surface wind, or from the 
1000 mb surface being a fictitious atmosphericlevel in many cases. 
Another explanation is that non-conservative influences on the 
temperature, Such as surface heating and turbulent mixing, are most 


intense below 850 mb. [2] 
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This particular evaluation method considered grid points in a 
fixed coordinate system, an Eularian approach. A Lagrangian 
approach ought to give better results, This was tested by selecting 
"identifiable points" midway between the baroclinic-zone 
boundaries and checking the time variations of YO against FG at 
these points. The results for the area in Figure 33, are given in 
Table 4. The most significant change from the results of the 
previous method is the increased verification for positive and 
negative FG at 850 and 700 mb. 

When the relation of FG to frontal movement was checked for 
all levels and for both temperatures fields, no correlation was 
observed. Centers of FG naturally moved in accordance with 
fronts, but gave no indication of future movement. 

In summary, a review of the two tables shows that FG 
(especially negative FG) is best verified at 850 and 700 mb. In order 
to improve the relation of frontogenesis to changes in Y6@, it was 
concluded that the non-conservative effects of the free atmosphere 
must be taken into consideration. 

8. SHVF Analyses 

The shear of the geostrophic wind component parallel to 
isentropes on a constant pressure surface (SHVF) (see equation 4) 
was obtained in units of 0.1 meters/sec/(100km). Figure 34 
schematically shows the vector geostrophic wind component 


parallel to isentropes (VF), as indicated by arrows, in four shear 
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producing situations. As with the other wind parameters, SHVF 
was computed at grid points for the five mandatory levels. 

Figure 36 shows the result over the North American area at 700 mb. 
Values of +8.0 meters/sec or approximately 16 knots/(100km) were 
the extremes observed, as found in the vicinity of GG®@ fronts at 
300 mb. The peak values decreased to +2.5 at the 1000-mb level, 
partially owing to decreased values of Vg. 

In order to determine the relation between SHVF and GG®@ 
fronts and baroclinic zones for the entire 63 x 63 grid, the following 
relationships were checked for both the processed and observed 
temperatures products. First, grid points within the baroclinic 
zone were categorized as to positive (cyclonic) or negative (anti- 
cyclonic) SHVF; secondly, all GG@ fronts were checked for 
algebraic sign of SHVF. Similarly, the USWB fronts were also 
checked at the only available constant pressure level, 850 mb. 

Table 5 shows the results of the first method for 00Z and 12Z 
20 August 1965. Only a slight overall tendency for grid points to 
experience positive SHVF rather than negative SHVF is noted. No 
large variation from this insignificant positive tendency exists 
at any level for either map time or with either temperature field. 

The results of the second method are shown in Table 6. Itis 
evident that there is a distinct tendency for fronts to lie mostly 
in positive SHVF at lower levels and predominately in negative 


SHVF at higher levels. This tendency shift from positive to 
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negative SHVF with height was observed at both times and with 
both temperatures, with only minor exceptions. Although it was 
not quantitatively investigated, it was noted qualitatively that 

the rear boundary of the baroclinic zone was characterized by 
positive SHVF at all levels. Table 7 shows an approximately 

equal distribution of USWB 850emb fronts in areas of positive and 
negative SHVF, with a larger variation for the observed-temperature 
data. This distribution compares well with GG@ fronts. 

The VF field was observed to be a combination of situations 
one and two of Figure 34, That is, maximum VF within the 
baroclinic zone produces negative SHVF at the warm-air boundary 
of baroclinic zones and positive SHVF at the cold-air boundary 
of baroclinic zones. VF maximum tends to shift in position from the 
vicinity of 95 to 95 with increasing height. Figure 35, showing 
VF plus a line of maximum VF and the baroclinic zones, and 
Figure 36, showing SHVF with respect to the GG@ baroclinic zones, 
illustrate this situation. 

According to the common school of thought for fronts, the 
baroclinic zone is characterized by large positive shear. The 
SHVF patterns and the GG@ baroclinic zones miss this characteri- 
zation by a small distance (less than one mesh length) at lower 
levels, but by a large distance (one half the width of the 
baroclinic zone) at higher levels. Maximum values of VF may be 


mislocated by as much as one half mesh length, due to the ability 
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of locating maximum values only at grid points, when in certain 
situations they should be more properly located between grid 
points. Although this lack of precision in locating VF maxima 
may have some effect on the SHVF-GG®@ baroclinic zone relation- 
ship at lower levels, it is too small to have a significant effect 
at higher levels, where the baroclinic zone has a width of two to 
four mesh lengths. So, it is concluded that here negative SHVF is 
characteristic of upper-level baroclinic zones and fronts, while 
positive SHVF dominates low-level baroclinic zones and fronts. 

To summarize, it is noted that the line of maximum VF within 
the baroclinic zone determines the peculiarities of the SHVF- 
baroclinic zone relationship for all levels tested. Also, the 
vertical variation of shear may be further evidence of dissimilar 
structure in fronts at low and high levels. 
9. Final Conclusions and Recommendations 

Except for the 1000-mb level, high latitudes, and the Asiatic 
continent, the processed temperature is considered to be a good 
approximation to the observed temperature, when moisture effects 
are considered. Since observed en verstices Suppress intensity 
of fronts over water surfaces at lower levels and is a poor front 
indicator at 1000 mb, it is considered inferior to GG@p at lower 
levels. At the three upper levels (700, 500, and 300 mb), GG@p 
and GG@8o are comparable. Perhaps a test using observed virtual 
temperature is now in order, as this would allow a better evaluation 


of the FNWF virtual temperatures products. 
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Of the three wind parameters studied, VNF manifests itself 
as being the most promising, especially at 850 mb,using processed 
temperatures. Its ability to indicate instantaneous front movement 
and its time consistency is suggestive of incorporating VNF into 
the numerical frontal analysis systems. 

Due to the numerous non-conservative effects present in the 
atmosphere, of which none were taken into account in determining 
frontogenesis, FG was found to be of limited value as an indicator 
of changes in intensity of baroclinic zones. 

Before stating positively that fronts lie predominately, although 
not: completely, in positive SHVF at lower levels and primarily in 
negative SHVF at higher levels, it is recommended that SHVF be 
recomputed using smaller grid lengths and/or utilizing a smaller 
computation interval. 

It would be improper to assume that the results of this study 
are applicable in general, since only two observation times of a 
Summer date were used. Therefore a similar study for a winter 
date would be beneficial in establishing the generality of the 
results obtained. 
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TABLE 1 


FNWF PROCESSED TEMP MINUS OBSERVED TEMP: 002Z 20 AUG 1965 
0.0 to 0.0 to +1.0to ~-1.0 to +3.0 to -3.0 to >+6.0 <-6.0 


+0.9 =(o +2.9 =, 48, +5.9 -5.9 
300 mb 
South 30 ie 45 6 5 1 
Central 11. 24 26 16 l 2 
North 22 25 20 25 Z 6 
500 mb 
South 47 ZS 16 8 3 1 
Central aly Sie 15 6 ] ] 
North 36 OW), 18 11 3 
700 mb 
South 43 14 4] 2 
Central 29 16 46 Saco Z i Oro 
North 29 22 Sil 2 51449 aD 
850 mb 
South 17 12 54 6 9.5 l wu ..5 
Central 20 i3 47 8.5 8.5 5) 
North 29 22 33 1] 3 Z 
1000 mb 
South 14 7 42 4 23 i 9 
Central Bh 1] 23 9 23 4 17 2 
North 17 13 26 AZ Leo 4.5 9 
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TABLE 2 


FNWF PROCESSED TEMP MINUS OBSERVED TEMP: 00Z 20 AUG 1965 
0.0 to 0.0to +1.0to -1.0to +3.0to -3.0 to >+6.0 €-6.0 


tO i 0.9 geet 2.9 ie +5.9 =o 
300 mb 
N.America 22 26 36 eZ 1 ZAS 
Pacific 38 23 28 8 Zn. 0.5 
Asia 24 20 27 17 if = 
Europe 28 19 39 9 4 1 
Atlantic 29 18 Sis) 14 3 
200 mb 
N.America 43 27 Zi 9 
Pacific Se: 27 14 4.5 ees 
Asia oul 32 de. 18 2 
Europe 40 26 21 6 7 
Atlantic 44 38 ie 7 
700 mb 
N.America 34 ln3 47 S S 
Pacific 38 24 30 6 1 1 
Asia 28 LL 48 6 4 OF OS 
Europe 37 21 ou 9 2 1 
Atlantic 39 6 49 Z 4 
850 mb 
N.America 18 LZ 54 8 Ogee Lato 
Pacific 19 6 63 3 S25 Og, 
Asia 16 19 S15, 18 7. 4.5 
Europe on eae | 30 Ta oa ero 
Atlantic 14 8 64 11 Lao tO 
1000 mb 
N.America 14 8 34 6 26 1 AT} 
Pacific ZS 14 41 9 eae 6 es iL, 1 
Asia 6 27 ZA 7 24 5 22 3 
Europe Ja 8 26 4.5 ol eS 17 
Atlantic zZ az 39 8 26 1 2 
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TABLE 5 


SHVF 
00002 12002 
POSITIVE NEGATIVE POSITIVE NEGATIVE 
(in %) (in %) (in %) (in %) 
PROCESSED TEMPERATURE 
1000 mb oy 43 58 42 
850 mb oy 43 54 46 
700 mb ov 43 46 54 
500 mb So 4S 49.9 SUeS 
300 mb S) 49 48 92 
OBSERVED TEMPERATURE 
1000 mb 48 OZ 
850 mb 99 4S a2 48 
700 mb a2 48 
500 mb 56 44 44 56 
300 mb 49 ot 
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TABLE 6 


SHVF 
0000Z 1200Z 
POSITIVE NEGATIVE POSITIVE NEGATIVE 
(in %) (in %) (in %) (in %) 
PROCESSED TEMPERATURE 
1000 mb iS, 21 60 40 
850 mb ac 48 04 46 
700 mb She 69 22 78 
900 mb 21 fs) Le 88 
300 mb 14 86 02 98 
OBSERVED TEMPERATURE 
1000 mb 98 42 
850 mb as 4l 52 48 
700 mb Zo 74 
900 mb lef 83 15 85 
300 mb 07 gs 
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TABLE 7 


SoUsmb aot 


PROCESSED TEMPERATURE 


POSITIVE NEGATIVE 
(in %) (in %) 
00Z 50 50 
27, 58 42 


4S 


OBSERVED TEMPERATURE 


POSITIVE 
(in %) 


40 


62 


NEGATIVE 
(in %) 


60 
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